arXiv: 1505.04762vl [math.PR] 18 May 2015 


Expected number of real zeros for random Freud orthogonal 

polynomials 

Igor E. Pritsker®'’*, Xiaoju Xie^ 

^‘Department of Mathematics, Oklahoma State University, Stillwater, OK 74078, USA 


Abstract 

We study the expected number of real zeros for random linear combinations of orthogonal 
polynomials. It is well known that Kac polynomials, spanned by monomials with i.i.d. 
Gaussian coefficients, have only (2/7r + o(l)) logn expected real zeros in terms of the degree 
n. On the other hand, if the basis is given by orthonormal polynomials associated to a hnite 
Borel measure with compact support on the real line, then random linear combinations have 
n/y/3 + o{ n) expected real zeros under mild conditions. We prove that the latter asymptotic 
relation holds for all random orthogonal polynomials on the real line associated with Freud 
weights, and give local results on the expected number of real zeros. We also show that the 
counting measures of properly scaled zeros of random Freud polynomials converge weakly to 
the Ullman distribution. 

Keywords: Polynomials, random coefficients, expected number of real zeros, random 
orthogonal polynomials, Freud weights. 


1. Background 


Problems on the number of real zeros for polynomials with random coefficients date 
back to 1930s, and they are considered as some of the most classical in the area of random 
polynomials. These original contributions dealt with the expected number of real zeros 
E[i\y (R)] for polynomials of the form Pn{x) = Efc=o where {ck}^^Q are independent and 
identically distributed random variables. Apparently the hrst paper that initiated the study 
is due to Bloch and Polya [^, who gave an upper bound E[A„(R)] = 0{y/n) for polynomials 
with coefficients selected from the set {—1,0,1} with equal probabilities. Further results 
generalizing and impro ving that estimate were obtained by Littlewood and Offord 


|2fil |. Erdos and Offord 
asymptotic result 


lOj and others. In particular, Kac [l8| established the important 


E[A„(R)] = (2/7r + o(l)) logn as u —)■ cxo. 


for polynomials with independent real Gaussian coefficients. More precise forms of this 


asymptotic were obtained by many authors, including Kac [l9j, Wang [3^, Edelman and 
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Kostlan 
of Wilkins 


It appears that the sharpest known version is given by the asymptotic series 


35| . Many additional references and further directions of work on the expected 


number of real zeros may be found in the books of Bharucha-Reid and Sambandham [l|, 
and of Farahmand 0. In fact, Kac [l^-jl^ found the exact formula for E[iV„(M)] in the 
case of standard real independent Gaussian coefficients: 


where 


E ^ ^ W \ 

TT Jo A{X) 


A{x) = B{x) = ''^jx'^^ ^ and C{x) = ^ j 




j=0 


i=i 


i=i 


In the subsequent paper Kac [2^, the asymptotic result for the number of real zeros was 
extended to the case of uniformly distributed coefficients on [—1,1], Erdos and Offord 10 
generalized the Kac asymptotic to Bernoulli distribution (uniform on { —1,1|), while Stevens 


32[ considered a wide class of distributions. Finally, Ibragimov and Maslova extended 


the result to all mean-zero distributions in the domain of attraction of the normal law. 

We state a result on the number of real zeros for the random linear combinations of rather 
general functions. It originated in the papers of Kac 18|- 20|, who used the monomial basis. 


and was extended to trigonometric polynomials and other bases, see Farahmand 0 and 
Das 0 - 0 . We are particularly interested in the bases of orthonormal polynomials, which 
is the case considered by Das j^. For any set E C C, we use the notation Nn{E) for the 
number of zeros of random functions fll.ll) (or random orthogonal polynomials of degree at 
most n) located in E. The expected number of zeros in E is denoted by E[iV„(-S)], with 
E[iV„(a, b)] being the expected number of zeros in (a, b) C M. 

Proposition 1.1. Let [a,b] C M, and consider real valued functions gj{x) G C^{[a,h\)^ j = 
0,... ,n, with go{x) being a nonzero constant. Define the random function 


Gn{x) = ^Cjgfix), 
j=0 


( 1 . 1 ) 


where the coefficients Cj are i.i.d. random variables with Gaussian distribution a > 

0. If there zs M G N such that G'.^{x) has at most M zeros in {a,b) for all choices of 
coefficients, then the expected number of real zeros of Gn{x) in the interval {a,b) is given by 


E[iV„(a,6)] 


1 ^A{x)C{x) - Bfix) 

TT X A{x) 


where 


( 1 . 2 ) 


n n n 

Mx) = '^g^{x), B{x) = '^gj{x)gj{x) and C{x) = '^[gj{x)]‘^. (1.3) 

j=0 j=l j=l 
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Clearly, the original formula of Kac follows from this proposition for gj{x) = x\ j = 


0,1,... ,n. For a sketch of the proof of Proposition 11.11 see |27|. We note that multiple 


zeros are counted only once by the standard convention in all of the above results on real 
zeros. However, the probability of having a multiple zero for a polynomial with Gaussian 
coefficients is equal to 0, so that we have the same result on the expected number of zeros 
regardless whether they are counted with or without multiplicities. 


2. Random orthogonal polynomials 

In this paper, we consider the Freud weights 

W{x) = X e M, 

where c > 0 and A > 1 are constants [l^. For n > 0, let 

Pn ( 2 :) = Pn x) = 7„x” + .., 

denote the nth orthonormal polynomial with 7 „ > 0, so that 


PnPmW^ = K 


Using the orthonormal polynomials {pj}’^o as the basis, we consider the ensemble of random 
polynomials of the form 


Pnix) = 


( 2 . 1 ) 


-\x) = ^CjPj{x), n e N, 

j=0 

where the coefficients cq, ci,..., c„ are i.i.d. random variables. Such a family is often called 
random orthogonal polynomials. If the coefficients have Gaussian distribution, one can 
apply Proposition 11.11 to study the expected number of real zeros of random orthogonal 
polynomials. In particular. Das considered random Legendre polynomials, and found 
that E[A^„(—1,1)] is asymptotically equal to nj\f%. Wilkins improved the error term 
in this asymptotic relation by showing that EWW—l, 1)] = nj^Z + o(n^) for any £ > 0. 
Related results were obtained by Farahmand jl2l. Il3|. For random Jacobi polynomials. Das 
and Bhatt [sj concluded that 1,1)] is asymptotically equal to n/\/3 too. They also 

stated estimates for the expected number of real zeros of random Hermite and Laguerre 
polynomials, but those arguments contain signihcant gaps. The results of this paper provide 
detailed information on the expected number of real zeros for random polynomials spanned 
by the Freud orthogonal polynomials. In particular, they cover the case of random Hermite 
polynomials. Lubinsky and the authors also recently showed 27| that if the basis is given by 
orthonormal polynomials associated with a hnite Borel measure compactly supported on the 
real line, then random linear combinations have n/\/3 + o(n) expected real zeros under some 
mild conditions on the weight. Interesting computations and pictures of zero distributions of 
random orthogonal polynomials may be found on the chebfun web page of Trefethen [33 
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For the orthonormal polynomials {pj{x)}’^Q, dehne the reproducing kernel by 


n—1 


Kn{x,y) = 

i=o 


and the differentiated kernels by 


n—1 


= '^pf\x)p^l\y), k,l eNU{0}. 

j=0 

The strategy is to apply Proposition 11.11 with gj = pj, so that 

A{x) = Kn+i{x,x), B{x) = K^_l^^{x,x) and C{x) = K^^^{x,x). (2.2) 

We use universality limits for the reproducing kernels of orthogonal polynomials (see Levin 
and Lubinsky [^-[^), and asymptotic results on zeros of random polynomials (cf. Pritsker 
3) to give asymptotics for the expected number of real zeros for a class of random orthog¬ 


onal polynomials associated with Freud weights. 

Theorem 2.1. Let W{x) = he a Freud weight on R, where c > 0 and A > 1 are 

constants. Then the expected number of real zeros of random orthogonal polynomials fl2.ip 
with Gaussian coefficients satisfy 


lim -E[iV„(R)] = 

n^oo n 


(2.3) 


We note that asymptotic relation fl2.3p is new even in the classical case of Hermite weight 

1 2 1^—I 

W{x) = e“ 2 ^ . Theorem 12.11 is a combination of two results on zeros of random orthogonal 
polynomials given below. Define the constant 


7a = 


r(j)r(f: 

2r(i±i) 


and the contracted version of 


P’(s) := P„ia„s), IJ 6 N, 


(2.4) 


(2.5) 


where On = 'yfc~^n^ is a positive number. 

For any set P C C, we use the notation N*{E) for the number of zeros of random 
functions Pf{s) located in E. The expected number of zeros of Pf{s) in E is denoted by 
E[iV*(P)], with E[iV* ([a, 6])] being the expected number of zeros in [a,b] C R. 

Theorem 2.2. Let Wfx) = he a Freud weight on R, where c > 0 and A > 1 are 

constants. If[a,b] C (—1,1) is any closed interval, then 


lim -E[iV* ([a,6])] = ^p^([a,6]), 

n->.oo n y 3 


( 2 . 6 ) 
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where the measure is given by 





ds, 


■S ^ [“!) !]• 


31 


Note that fiw is the weighted equilibrium measure for the weight w{x) = e 


— p-7Ahl'^ 


on. 


see 


and the next section for details. This measure is often called the Ullman distribution. 
Dehne the normalized zero counting measure for the scaled polynomial 

F*(s) of fl2.5p . where {zk}1=i are its zeros, and 5z denotes the unit point mass at z. We can 
determine the weak limit of for random polynomials with quite general random coefficients 


Theorem 2.3. If the coefficients {cfcj^o fcmdom orthogonal polynomials fl2.ll) are com¬ 
plex i.i.d. random variables such that E[|log|co||] < oo, then the normalized zero counting 
measures Tn for the scaled polynomials Pf{s) converge weakly to /i^ with probability one. 


Closely related results on the asymptotic zeros distribution of random orthogonal poly¬ 
nomials with varying weights were proved by Bloom and Bloom and Levenberg j^, but 
they are not directly applicable to our case because of different normalization. Theorem 12.31 
allows to find asymptotics for the expected number of zeros in various sets. In particular, 
we need the following corollary for the proof of Theorem 12.11 

Corollary 2.4. Suppose that the coefficients {cfcj^Q of random orthogonal polynomials fl2.ip 
are complex i.i.d. random variables such that E[|log|co||] < oo. If E G C is any compact 
set satisfying fi^idE) = 0, then 

Mm = (2.7) 

n—>-oo 77, 

where Nf{E) is the number of real zeros of Pf{s) in E. 


It is of interest to develop similar results for a more general class of weights W on the 
whole real line, defining orthogonal polynomials pk. Another direction of further work is 
related to relaxing conditions on random coefficients c^, e.g., W considering probability 
distributions from the domain of attraction of normal law as in 16. 1^ . 


3. Proofs 

Our proofs require detailed knowledge of potential theory with external fields generated 
by Freud weights, see 22| and (3lj |. 

Let IF be a continuous nonnegative weight function on R such that IF is not identically 
zero and hm| 2 .|_,,oo |a;| lF(a;) = 0. Set Q{x) := — loglF(a;). The weighted equilibrium measure 
fiw of is the unique probability measure with compact support Sw = supp fiw C R that 
minimizes the energy functional 


ns] = 


\og\z — t\dfi{t)dp.{z) + 2 Q dpi 
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amongst all probability measures /i with support on R. It satisfies 


log 


and 


log 


z —1\ 
1 


z — t\ 


+ Q{z) — C*, Z E. S\y, 


dfiwit) + Q{.z) > C, 2 ; G R, 


where C is a constant. We consider the general Freud weights W{x) = e where c > 0 
and A > 1 are constants, and the normalized weight w{x) = where 7 a is defined by 

02.dp . In the latter case, the weighted equilibrium measure is given by 




^ “'hi 


,A-1 




■dy]d3, s e [-1,1], 


by Theorem 5.1 of [3ll, p. 240]. 

For a weight function W{x) = the Mhaskar-Rakhmanov-Saff numbers 


are defined for n > 1 by the relations 

1 r 


n = 


a_„ < 0 < m 


xQ'{x) 


TT Ja_„ - a-n){an - x) 


dx 


and 


0 = 


Q'{x) 


TT Ja_„ \/{x- a_n)(a„ - x) 


dx. 


Since Q{x) = c\x\ is even, we have a_r, = —Ur,. Existence and uniqueness of these numbers 
is established in the monographs j^, j^, 31], but goes back to earlier work of Mhaskar, 
Saff, and Rakhmanov. One illustration of their role is the Mhaskar-Saff identity: 


||FW|h 


I |-PhF| |Loo([-a„,a„])5 


which is valid for all polynomials P of degree at most n. It is known that the Mhaskar- 
Rakhmanov-Saff number associated with the Freud weight W(x) = is given by 

i -- - 

an = 7a c 


See page 308 of ^ for further details. We define the Mhaskar-Rakhmanov-Saff interval 
as := [—an,a„]. The linear transformation 

X 

Ln{x) = —, x G R, 
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maps An onto [—1,1]. Its inverse is 

= a„s, s e M. 

For e G (0,1), we let 

Jn{e) = + e,l - e] = (1 -£)[-a„,a„]. 

Then the equilibrium density is dehned as 

+ x){an - x) Q'{s) - Q'{x) ds 


an{x) = 


TT^ 


^ ^ \/ {o,n + s){an — s) 


, X e Ar, 


The equilibrium density satishes 22|, p. 41]: 


log -fCrn(s) ds + Q{x) = C, X e An, 


X — s 


and 


log -j-rO'n(s) ds + Q{x) > C, X E 


\x — s 


Note that the measure an{x)dx has total mass n: 

/ an 

crn(x) dx = n. 

an 

We also dehne the normalized version of an as follows: 


(T*(s) := —an{ans), s e [-1,1]. 
n 


Note that 


<(s) ds = 1. 


'-1 


For details on an, one should consult the book [2^ by Levin and Lubinsky. 

Lemma 3.1. For a Freud weight W{x) = , where c > 0 and A > 1 are constants, the 

normalized equilibrium density satisfies 


That is, 


(T*(s) ds = djjini^s) for all n Efi, s E [—1,1]. 

dy for alln E'N, s E [—1,1]. 


A r 


<(s) = - I ^ 
^ J\s\ 


2 _ g2 
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Proof. Recall that satisfies 


log y-- -7(Tn{y) dy + C |x|^ = Ci, X e [-On, ttn], 


and 


\x - y\ 


log 


lx - y\ 


(Tn{y) dy + C |x|'^ > Cl, X G M, 


^ 1 1 


where = 'yfc The changes of variables x = a^s and y = Ont rednce the above 

relations to 

1 .,A 


»1 


and 


Li |s-t| 
rl 1 


<(t) dt + 7 a |s| = Ca, s e [-1,1], 


/ log I-TT 

'-1 |s-^l 


(T*(t) dt + 7a |s|^ >C2, s G M. 


Invoking Theorem 3.1 of j3l|, P- 43], we dednce that 

(T*(s) ds = d+w{s) for all n G M, s G [—1,1]. 
On the other hand, recalling that 


dpio(s) I 


X /•! 2/A-I 



:dy\ds, s G [—1,1], 


Id 


we obtain that 


<(s) = 


A C y^-^ 


^ J\s\ \/y‘^ - 


\/y‘^ - 

dy for all n G N, s G [—1,1]. 


□ 


Proof of Theorem \2.S[ In this case, hT(x) = G M, where c > 0 and A > 1 are 

constants. The strategy is to apply Theorem 1.6 of 23(|. It states that for all r, s > 0 and 
any e G (0,1), we have 


W‘^{x)Kif’"\x,x) 


E 

j=0 



{an{x)Y+^+^ 
nniformly for x G JnY), where 


k=0 



E U V.- 


.j+k 


Q'(x) 


(Xn X 


r-\-s—j—k 


+ o(l) as n —)■ cxo, 

(3.1) 


P,k 


0, j + k odd, 

(_l)0-0/2_^^ j + A: even. 


That is, nniformly in x G Jn+iY), 

W\x)K+l (x,x) 

(yn+l{x) 


= 1 + o(l) as n ^ oo. 


















and 


WHx)K;°i\Ux,x) _ Q'(x) 


(an+l(x))2 

dl>l) 


0-n+l(x', 


+ o(l) as n oo, 


W^(x)K!;['(x,x) ( Q'(x) >r^ 

Applying fll.2p with fl2.2p . we obtain that 


A'i+V {x, x) 




dx 


(3.2) 


for any closed interval [/, q] C Jn+i{s) (/, q may depend on n). Now, nniformly for x G Jn+i{e), 


and 


Kj^+i {x, x) 
kI^+i {x, x) 

ATj+y (x, x) 

^A:i+¥ {x,x)^ 






= ‘^^+.W|-^ + o(l)) (1+0(1))-" 


= (^n+AX) 


(Tn+l{x) 

Q\x) 


a^+iix) 




We thns obtain that 

^E[W([/,g])] 


'x Ji n 


\ 




{Ml) +f+ 0(1) , , SAAi) +^0(1) + 0(1) 


1 + 0 ( 1 ) 


-o5+i(a:) 


1 + 0 ( 1 ) 


dx 


i + f ^)h(i)+^ 0 ( 1 ^ 7 ^* 

TV Ji n y 3 Vc^n+i(ao)/ cr„+i(x) 


as n —)■ cx). 


Note that the nnmber Nn{E) of real zeros of Pn{x) in E eqnals the nnmber N*{E*) of 
real zeros of P^(s) in E* := E/an+i, since A„+i is a bijection. Since [a,b] C (—1,1) is 
a closed interval, we have that [a„+ia, a„+i6] C Jn+i{£) = (1 — £()[—a„+i, a„+i] provided 
max{|a|, |6|} < 1 — e for some constant e G (0,1). Hence 


^E[N*{[a,b])] = ^E[Nn{[an+ia, ttn+ib])] 


(3.3) 


1 

TT 




' ^n + l<^ 

n + 1 /■' 


—(/t + 

n y 3 Vcrn+i(a;)/ o-n+i(a;) 


vm 


<+i(s) 




TT^ 






++ 1 1 o(l) ds. 
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We show that |Q'(x)/ct„+i(x)| < C on x G [a„+ia, a„+i6] for some constant C > 0 indepen¬ 
dent of n. Recall that [a„+ia, a„+i6] C Jn+i{s) = (1 — ^)[— a„+i], for some £ G (0,1). 
It is clear that Q'{x)/an+i{x) is an odd function of x G [a„+i, 0) U (0, a„+i], so that we only- 
need to consider the interval (0, a„+i(l — e)]. First note that for x G (0, a„+i), 


Q\x) 

Crn+l{x) 


^A-1 


l-l/A 1/A X 

= " " (n+l)i-VA 



y 


A-l 


-1 


— x^/a: 


2 

n+l 


dy 


Since A > 1, we can estimate the integral in the above formula for x G (0,an+i(l 
follows: 


y 


A-l 


dy > 


y 


A-l 


Aa;/a„+i \/^ Jx/a^+i \/W 

Using this estimate, we see that for x G (0, a„+i(l — e)], 


dy = 


1 - (x/an+i)^ ^ 

A-l 


£)] as 


Q\x) 

CTn+lix) 

< „1-1/A 1/A (Qrt+l(l - g))^ ^ A-l 

(n + l)i-iA i_(i_£)a-i 

^ 7A7r(A-l)(l-£)^-^ ^ 

l-(l-£)^-l 




X 


A-l 


A-l 

T'a ^(^+l)l-l/Al_(^/a^^^)A-l 


Applying (I3.3p and Lemma 13.11 we obtain that 


1 

n 


■E[iV:([a,6])] 



1 + 0(1) 

a/3 

1 + 0 ( 1 ) 

a/3 



<+i(s) ds 



djJ^w ('^) • 



+ o(l) ds 


To complete the proof, we pass to the limit as n ^ oo: 

1 1 /■^ 1 
lim -E [A"* ([a, b])] = -^ dfx^{s) = ([a, b]). 

n^oo n y^Ja a/3 

□ 

Proof of Theorem A2.d\. Following j^, we call a sequence of monic polynomials {Qn}'^=i, 
with deg(Q„) = n, asymptotically extremal with respect to the weight w if it satishes 


lim W'w'^Qn 

n—^oo 


II 1/n 

Hr 


= e 


—Fw 


1 


where || • ||r is the supremum norm on M and = log 2-|- 1/A is the modihed Robin constant 


corresponding to w, see [Ml, P- 240]. Theorem 4.2 of [3l|, P- 170] states that any sequence 
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of such asymptotically extremal monic polynomials have their zeros distributed according to 
the measure Namely, the normalized zero counting measures of Qn converge weakly to 
/i^. We show that the monic polynomials 

Qni^) ■= Pn{x)/{Cnlnal), U G M, 


are asymptotically extremal in this sense with probability one, so that the result of Theorem 
12.31 follows. 

Using orthogonality, we obtain for polynomials dehned in fl2.ip that 


\Pn{x)\^W‘^{x) dx = '^ |Cfcp. 


Hence 


max Icfcl < 

0 <fc<n 


fc =0 


1/2 


\Pn{x)fW^(x) dx ] < (n + 1) max Icfcl. 

' 0<fc<n 


Lemma 4.2 of [30|] (see (4.6) there) implies that 


lim 

n—>-oo 


\ l/(2n) / \ 1/n 

\Pnix)\‘^W‘^(x) dx] = lim f max Icfcl 1 =1 

J n^oo yo<k<n 


with probability one. Applying the Nikolskii-type inequalities of Theorem 6.1 and Theorem 
6.4 from j^, we obtain that the same holds for the supremum norm: 


lim 


= lim \\PnW\\P^ = 1 

n—>-oo 


with probability one. Recall that the leading coefficients of orthonormal polynomials Pr^ 
satisfy 

hm7y"ni/" = 2c'/^;'/V/" 


by Theorem 1.2 of 31^. 362]. We also use below that hm„_,.oo = 1 with probability 

one by Lemma 4.2 of 3^. It follows that 


lim = lim \\w'^P*\\P'' lim |cn7n| |an| lim 

n—>-00 n—>-00 n—>oo n—>-00 \ / 


-1 


-1 


7-l/A^l/A^_l/A7l/A^| ^ g-(log 2 +l/A) ^ g-F™ 


□ 


Proof of Corollary \2.4\ Consider the normalized zero counting measure Tn = ^ J2k=i^zk 
the scaled polynomial Pf{s) of fl2.4p . where are the zeros of that polynomial, and 

5z denotes the unit point mass at Theorem 12.31 implies that measures converge weakly 
to pu, with probability one. Since pwidE) = 0, we obtain that Tn\E converges weakly to 
Pw\e with probability one by Theorem 0.5' of [2l| and Theorem 2.1 of j^. In particular. 
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we have that the random variables Tn{E) converge to /i^(-E) with probability one. Hence 
this convergence holds in sense by the Dominated Convergence Theorem, as Tn{E) are 
uniformly bounded by 1, see Chapter 5 of [l^. It follows that 

lim E[\Tn{E) - fi^{E)\] = 0 

n^oo 

for any compact set E such that fiw{dE) = 0, and 

\E[Tn{E) - iJ.w{E)]\ < E[\Tn{E) - iJ.nj{E)\] -^0 as n -)■ oo. 

But E[r„(T^)] = E[N*{E)]/n and E[fi^{E)] = fi^{E), which immediately gives fl2.7p . □ 

Proof of Theorem \2.1[ Theorem 12.21 gives that 

lim -E [iV* ([a, b])] = ^/i^([a, b]) 
n^QO n i/3 

for any interval [a, b] C (—1,1). Note that both E [iV* {H)] and ^w{H) are additive functions 
of the set H. Moreover, they both vanish when H is & single point by (12.7p and the absolute 
continuity of with respect to Lebesgue measure on Sw = [—1,1]. Hence (12.71) gives that 

lim -E [K (M \ (-1,1))] (-1,1)) = 0. 

n—>-oo fi 

It now follows that 

lim -E[iV*(R)] = 

n^oon y/3 y/3 

To complete the proof, observe that Nf(R) = so that E[iV*(M)] = E[A^„(M)], since 

Ln+i is a bijection for each hxed n. Therefore fl2.3p is proved. □ 
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